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ABSTRACT 

New, more detailed calculations of the emission spectra of the He-like Ka complex of 
calcium, iron and nickel have been carried out using data from both distorted- wave and 
R-matrix calculations. The value of the GD ratio (an extended definition of the G ratio 
that accounts for the effect of resolved and unresolved satellite lines) is significantly 
enhanced at temperatures below the temperature of He-like maximum abundance. Fur- 
thermore it is shown that satellite lines are important contributors to the GD ratio 
such that GD /G > 1 at temperatures well above the temperature of maximum abun- 
dance. These new calculations demonstrate, with an improved treatment of the KLn 
(n ^ 3) satellite lines, that Ka satellite lines need to be included in models of He-like 
spectra even at relatively high temperatures. The excellent agreement between spectra 
and line ratios calculated from R-matrix and distorted-wave data also confirms the 
validity of models based on distorted-wave data for highly charged systems, provided 
the effect of resonances are taken into account as independent processes. 

Key words: atomic data, atomic processes, line: formation, line: profiles, X-rays: 
general 



1 INTRODUCTION 



As d escribed in the literature (e.g. iGabriel fc Jordan! 
Il969al h the Ka emission of helium-like ions takes 
place via four lines, with the designations w, x, y 
and z: ls2p *Pi -» Is 2 % (w), ls2p 3 P 2 -> Is 2 1 S (x), 
ls2p 3 Pi -> Is 2 l So (y), ls2s 3 Si -> Is 2 ^0 (z). From these 
four lines, line ratios have been investigated for diagnostic 
purposes. One of these ratios, G, is sensitive to temperature 
and is defined as 

/(x) + /(y) + /(z) 



G = 



Kw) 



(1) 



where /(w) is the intensity of the w line in units of num- 
ber of photons per unit volume per unit time. For heav- 
ier elements, additional lines arising from transitions of 
the type ls21nl' — > ls 2 nl', where the upper ls21nl' states 
are autoio nizing, tend to complicate this simple spec- 
trum (e.g. lEdlen fc Tvrenl Il939l; IGabriel fc Jordan! Il969bl : 
iGabriell 1 19721 ; iMewe fc Schriiverl Il97ct ). The KLL satellite 



lines (which arise from configurations o f the type ls2 121') 
are designated with the letters a-v (see IGabriell 1 19721: th e 
most recent treatment is given bv lNahar fc Pradhanll200cj) . 
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Higher satellite lines arising from ls21nl', n ^ 3, are usu- 
ally not given separate designations. Often, astrophysical 
spectra can not be measured such that these satellite lines 
are adequately resolved; the result is what appears to be a 
broadened and redshifted w line according to the intensi- 
ties of the KLL lines in t o to within the Kg complex (e.g. 
Oelgoetz fc Pradhadl200ll: iHellier fc MukailliooilXu et all 



200' 



ISwartz fc Sulkanenl (|l993l ) first proposed a method of 
analysing emission spectra if the resolution was sufficient 
to resolve the spectra into two ranges, one corresponding 
to an energy range around the w line, the other including 
everything redward of this range. The range about the w 
line would include not only the w line itself, but also most 
of the satellite lines arising from the configurations ls21nl', 
n ^ 3. The redward range would include the bulk of the KLL 
satellite lines, in addition to the x, y and z lines. Thus, they 
proposed redefining G (referred to as Gs&s below) by taking 
the integral of the flux redward of some specified boundary 
line and dividing by the integral of the flux blueward of that 
same line. 

iBautista fc Kallmanl (|2000h considered the same ef- 
fect by including the intensity of the satellite lines as 
par t of the numerator in their calculation of G. Soon af- 
ter, [Oelgoetz^&^£radha3 (|200ll ) proposed a new ratio, GD, 
which included all the KLL satellite lines in the numerator 
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and all the satellite lines arising from higher shells in the de- 
nominator. For some elements, a weak KLL line is present in 
the area one would associate with the w line; additionally, 
for most heavier elements, some higher lines (which arise 
from configurations such as ls21nl' , n 3) have low enough 
energies such that they s hould be included with th e x, y and 
z lines in the numerator (Bclv -Dubau et al ] |l979al lbh. To im- 
prove on these earlier efforts, the GD line ratio is redefined 
in the current work as 



GD 



E H*) 

nin <E s <E b 



E 

E b <E s <E Ka 



I( S ) 



(2) 



where E b is the boundary line between the two energy 
ranges, Eicct-mm and -Eifa-max denote the energy range of 
the Ka complex, and E„ is the energy of a particular line, s. 
Thus, each sum includes the intensity of each line which has 
its centroid in the appropriate range. In the low tempera- 
ture limit, Gs&s w GD, but at high temperatures, Doppler 
broadening will cause the wings of lines near Eb to appear 
in the other range when computing Gs&s- 

In addition to being a temperature sensitive diagnostic, 
the G and GD rati os are also sensi t ive to the ionization 
state of th e plasma (|Pradhanl Il985bl ; lOelgoetz fe Pradhanl 
2001, 2004). While plasmas out of coronal equilibrium are 
not considered in this work, the results presented here have 
direct implications and utility to modeling those systems. 

Lastly, it should be noted tha t some recent work 
|Rana et all |2006l : iGirish et all 120071 ') omits the satellite 
lines from analysi s of Fe Kg observations on the basis of 
the argument that lOelgoetz fe Pradhanl (|200lh showed that 
the contributions from these lines can be neglected above 
the temperature of He-lik e maximum abundance. While 
lOelgoetz fc Pradhanl (120011) reported that GD « G in the 



range T e > 3.0 x 10 7 K (|Oelgoetz fe Pradhanll200ll , Fig. 2), 



they also showed that the satellites a re an important part 
of th e flux in this temperature range (|Oelgoetz fe Pradhanl 
l200ll . Fig. 1). In these earlier calculations, the contribution 
of satellite lines to the denominator and numerator of GD 
effectively cancelled each o ther out, resulting in GD ~ G. 
However, the calculations of lOelgoetz fe Pradhanl (|200ll ) did 
not treat the KLn (n > 3), satellite lines on par with the 
KLL lines. Specifically, the K LL lines were trea t ed rig orously 
according to the method of iGabriel fe Pagetl (|l972T l. while 
the KLn (n 3) lines were treated more approximately via a 
scalin g of ratios of autoioniza tion rates. Additionally, recent 
work (|Oelgoetz et al.l l2007al ) has shown that the cascade 
contr ibution to the recombination ra t es dMewe fe Schriiverl 
1 19781 ) used in lOelgoetz fe Pradhanl (|200ll ) diverges from 
the corresponding contribution calculated with more mod- 
ern distorted-wave and R-matrix methods at temperatures 
above the Fe He-like temperature of maximum abundance 
(T e >4xl0 7 K).For these reasons a new study which treats 
the KLn (n 3) satellite lines on par with the KLL satellite 
lines, and which is also based on more accurate atomic data, 
is warranted. 



2 THEORY 

The present work employed the General S pectral Model- 
ing (GSM) code jOelgoetzl |200d . see also lOelgoetz et~aH 



l2007al lbl). GSM is based o n the ground-stat e-only quasi- 
static approximation (e.g. iBates et al.l Il962l ). a common 
method for modelling low density plasmas such as found 
in astrophysics, which assumes that the ionisation balance 
portion of the model can be separated from a determination 
of excited-state populations. The rationale for this approxi- 
mation is two fold: first, the times scales for ionisation and 
recombination are much longer than the time scales for pro- 
cesses inside an ionisation stage, and second, the populations 
of the excited states have a negligible effect on ionisation and 
recombination. fSee lOelgoetz et al.ll2007bl . for a discussion 
of the validity of this approximation.) Thus, the first step in 
a GSM calculation is to solve the coupled set of ionisation 
balance equations given by 



dXi 
dt 



iv e (A' i+1 Q i+1 _ i (re)+x ; _ 1 c i _ 1 _ >i (r e )) 

+N e 2 (Xi +1 /3i +1 ^i(T e ) - jr,A-i_i(T e )) 
-XiN e (ai^i-!(T e ) + Ci^i +1 (T e )) , 



(3) 



where Xi is the total population in the i* 1 ionisation stage, 
N e the electron number density, T e is the electron tem- 
perature, C is a bulk collisional ionisation rate coefficient, 
P a bulk 3-body recombination rate coefficient, and a a 
bulk recombination rate coefficient (which includes radiative 
and dielectronic recombination). In general, photoionization 
and stimulated recombination are included as well, but as 
this work considers only collisional plasmas, the rate coef- 
ficients associated with these processes have been omitted 
from equation Q. 

Once the values of Xi have been determined, the 
ground-state-only quasi-static approximation then allows 
one to solve for the excited-state populations in a given ion- 
isation stage, with the approximation of treating the ioni- 
sation stages adjacent to the ionisation stage of interest as 
being entirely in the ground state. As the total population 
in the ionisation stage of interest and the two adjacent ion- 
isation stages are known, the excited-state populations can 
be determined by solving a modified version of the full set 
of collisional-radiative equations given by 



= JV„ ( E (Nt^ATe) - N^qf^n)) 



+ Xi+iCf+i t i->ij(T e ) — Ni,j E C ; °*^;-i,i(Te) 

i 

+ Xi-iaii!?. 1> i->i > j(T e ) - N Lj afj^i+i,j(T e )\ 

i(i>3) V 

^E^wP'.)) -Ni* E A f^ 

i J i(i<j) 

+ iV P ( E W,i£:f<T.) - N uq »-f(T e )) ) 



+ N a l J2( N ^f(Te)-N ljq fZf(Te)) 
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E D AI-off 



X, 



1,3 



(4) 
(5) 



where the variables are defined more or less as before; Nij 
is the population in the j state of the I ionisation stage, 
q cS is an electron-impact (de-)excitation effective rate coeffi- 
cient, g p ~ cff is a proton-impact (de-)excitation effective rate 
coefficient, g Q_cff is an alpha-particle-impact (de-)excitation 
effective rate coefficient, R Al ~ cB is an effective autoioniza- 
tion rate coefficient, and A eft is an effective radiative decay 
rate. Here we have used the "eff" superscript to denote the 
possible use of effective rate coefficients since GSM offers the 
option of treating some of the excited states as statistical 
conduits (using branching ratios) and others explicitly (e.g. 
lOelgoet dl2006l ; lOelgoetz et al.ll2007al fbh ■ Explicit states are 
those that appear in the set of coupled equations presented 
in equations Q and When all states within an ionisa- 
tion stage are treated explicitly, the "eff" superscript is not 
necessary since all of the rate coefficients represent direct 
processes only. When the statistical treatment is employed, 
the rate coefficients associated with the processes passing 
through statistical states are combined with the direct rate 
coefficients between explicit levels by summing over all the 
indirect paths through the statistical states. This process is 
simplified by the use of the collisionless transition matrix 
(CTM), T m ^j, which can be thought of as the probability 
that an ion in statistical state m will end up in an explicit 
state j, assuming that the time scale for collisions is very 
long when compared to the time scale for the spontaneous 
processes of autoionization and radiative decay. If Q repre- 
sents the set of explicit states, and i a state such that i Q, 
the CTM can be defined using the recursive expression 



E 



kgQ 
(E^E^Ej) 



Ai^i + 



+ 



(6) 



where Ti^k is the appropriate type of spontaneous rate (ei- 
ther radiative decay or autoionization) to connect states i 
and k . ft should be noted that if i £ Q the CTM is not 
meaningful, and as such is defined to be zero. The effective 
rate coefficient is then calculated by summing the direct 
rates, and the fraction (as determined by the CTM) of each 
indirect rate which contributes to an effective rate. For ex- 
ample, effective recombination (RR+DR) rate coefficients 
are calculated as 



eff Irr, \ RR- /Ti \ i T-iDC /T-i 



+ 



+ 



E 

(E m >E i 

E 

(E m >0) 
(E m >EJ 
(E m >E,-) 



(7) 



where D° c is a dielectronic capture rate and the sums take 
into account both radiative recombination and dielectronic 
capture followed by radiative cascade. It should be noted 
that there are terms in these sums that would be repre- 



sented by explicit resonances in R-matrix cross sections. 
Such an approach allows for the inclusion of resonances 
when perturbative (e.g. distorted-wave) cross sections are 
employed. This approach is sometimes referred to as the 
independent-process, isolated-resonance (IPIR) method (see 



Bates fc Dalgarnolll962l ; iGabriel fc Pagetll972l ; ICowanllf980| : 



garnollli: 

Badnell et all 19931 ) . In calculations that consider R-matrix 



data, care must be taken to exclude these terms from the 
summations in equation (J?) in order to avoid double count- 
ing the resonance contributions. 

Once excited-state populations have been calculated, 
the intensity of each line in the spectral region of interest is 
calculated according to 



I(l,j^l,k) = N ltj A^ k 



(8) 



Each line is then given a line shape corresponding to a ther- 
mal Doppler-broadened Gaussian profile. The total spec- 
trum (or emissivity), S, for a given photon energy, hv, can 
be expressed as 



S(hv) =hv^I s 



AE S V 2nkTi 



. g 2(AE s )' z kT i 



(9) 



where S is in units of energy per unit volume per unit time 
per energy interval, s ranges over the set of all included 
transitions in the desired energy range, AE 3 is the transition 
energy associated with a given line, and the ion temperature, 
Ti, is taken to be equal to the electron temperature. 



3 COMPUTATIONS 

As this work is concerned with steady-state plasmas, the 
solution to the coupled set of ionisa tion balance equations , 
Eq. Q, were taken to be those of iMazzotta et all (|f 9981 ) 
for all three elements (Ca, Fe, and Ni) considered in this 
work. Furthermore, as the cases considered are well within 
the low density limi t (N e = I0 10 cm -3 ), the approxima- 
tion IMazzotta et al.l l|f 9981 ) made in neglecting three-body 
recombination is valid. 

The present work considered multiple classes of mod- 
els for each of the three elements. Each model contains a 
different set of detailed atomic data. The first class, com- 
posed mostly of distorted- wave (DW) data (and denoted by 
Ni:DW, Fe:DW, and Ca:DW for the three elements), uses 
a set of data calculated entir ely by the Los Alamos suit e 
of atomic physics codes (e.g. lAbdallah et al.l [f994l . 1200 fl ). 
The CATS code was used to calculate the wave functions, 
energies, and dipole allowed radiative decay rates for all 
fine-structure levels arising from the configurations nl, lsnl, 
2ml', fs 2 nl, Is21nl', and Is31nl' with n < 10 and 1 < g, 
which span the H-like, He-like, and Li-like ionisation stages. 
The GIPPER code was used to calculate all autoionization 
rates and photoionization cross sections in the distorted- 
wave approximation, as well as collisional ionisation cross 
sections using a scaled-hydrogenic approximation which has 
been shown to agree well with distorted-wave results for 
highly charged systems. Distorted-wave cross sections for all 
electron-impact excitation transitions out of the lowest seven 
levels of the helium-like ionisation stage, as well as the Is 2 21 
complex of the Li-like ionisation stage were calculated with 
the ACE code. Cross sections for the remaining electron- 
impact excitation transitions were computed in the more ap- 
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proximate plane-wave Born approximation. Lastly, the non- 
dipole A values that give rise to the x and z lines, as well 
as a two- photon decay rate from ls2s 1 So — > Is 2 1 So used 
in obtaining the p opulations from equation s Q and ([5]), 
were obtained from iMewe fc Schriiverl (|l978l 'l . Proton- and 
alpha-particle-impact excitation rates between the He-lik e 
ls21 levels were also taken from IMewe fc Schriiverl (| 19781 ). 
The protons and alpha particles were taken to have the 
same temperature as the electrons, and to have densities 
of 0.77 and 0.115 times the electron density respectively 
(N p = 0.77iV e , N a = 0.115iV e ). The CATS level energies 
for the lowest seven levels of the He-like ionisation stage 
and the lowest three levels of the Li-like ionisation stage 
were replaced by values taken from t he NIST Atomic Spec- 
tra database (|Ralchenko et al]|2008l ). as were the energies 
for the KLL autoionizing levels for Li-like Ni and Fe. As the 
NIST database does not contain complete information for 
the autoionizing KLL levels of Ca, the level energies calcu- 
lated by CATS were retained for all Ca autoionizing states. 
All of the fine-structure levels arising from the Is, Is 2 , ls21, 
Is 2 21, and ls21nl' configurations with n ^ 10 and 1 ^ g were 
treated explicitly when solving for the excited-state popula- 
tions appearing in equations Q and ([5]). 

In the second class of models the electron-impact exci- 
tation, radiative decay, and both radiative and dielectronic 
recombination data in the DW model are replaced with data 
calculated using R-matrix (RM) methods, where such data 
are publicly avai l able. For Ni, the radiative decay rates of 
iNahar fc Eissnerl i|2006h and the unified recombination rates 
of iNaharl (|2005F were used to create the Ni:RM data set. As 
for Fe, two sets of R-matrix electron-impact excitation rates 
are available and are considered here. The first set, Fe:RM, 
in cludes th e electr on-impact excitation collision strengths 
oflPradhanl dl985al \. a su bset of the radiative decay rates of 
INahar fc Pradhanl (|1999h (where the initial state is a fine- 
structure level arising from the configurations lsnl where 
n ^ 4, 1 Sj / or 5 ^ n 10, 1 ^ p) and the correspond- 
ing su bset of the unified recombination rates of lNahar et al.l 
(|200ll ) (Is 2 Si/2 recombining into all fine-structure levels 
arising from lsnl where n 4, 1 / or 5 ^ n ^ 10, 1 ^ p). 
The second set, Fe:RM2, is identical to Fe:RM except that 
it uses the electron-imp act excitation collision strengths of 
IWhiteford et all (|200ll ). Lastly, one R-matrix type model 
is considered for Ca, Ca:RM, wh ich also incorpor ates the 
electron- impact excitation data of lPradhanl (|l985al ). 

The last class of models is an expansion of the second 
class by also incorporating autoionizati on rates calculate d 
from recom bination cros s sectio ns (e.g. INahar et al]|2007f ). 
Specifically, INahar et all l|2007h provided this type of data 
for Fe and Ni. These data have been combined with the 
Fe:RM and Ni:RM sets to make the Fe:RM+ and Ni:RM+ 
sets. The F e:RM2 data set, which i ncorporates the collision 
strengths of lWhiteford et all (|200ll ). has not been expanded 
into a Fe:RM2+ data set due to the good agreement (which 
is shown in the following section) between the Fe:RM and 
Fe:RM2 data set. As no data of this type are yet available 
for Ca, no model of this class is considered for Ca. 

In addition to constructing the models, the bound- 
ary line between the high energy and low energy section 
of each spectrum had to be chosen. As pointed out by 
ISwartz fc Sulkanenl (|l993l ) there is an energy gap that forms 
between the w line (and the satellite lines that blend with it) 



and the rest of the spectrum. This gap was found by inspec- 
tion, and the boundary energy, Eb, was chosen to be 3895 
eV, 6690 eV, and 7794 eV for Ca, Fe, and Ni respectively. 



4 RESULTS 

Figs. [TH3] display the calculated values of the G and GD 
ratios as a function of temperature for each of the models, 
along with plots of certain ratios that help to illustrate where 
the differences occur. 

Overall, the present calculations predict GD ratios 
that are significantly higher than the corresponding G 
ratios for all the models that are considered. This be- 
havi our is in qualitative agr e ement with previous stud- 
ies (ISwartz fc Sulkanenl 1 19931 ; iBautista fc Kallmanl |2000| ; 
lOelgoetz fc Pradhanl l200ll ); it should be noted that this 
more detailed study predicts a significantly greater value 
of GD below the temperature of maximum abundance than 
any of the previous studies. Additionally, the impact of satel- 
lite lines on the GD ratio keeps the GD/G ratio greater than 
one over a much bro a der ra nge than shown in the study of 
lOelgoetz fc Pradhanl (|200ll ). The principal reason for this 
behaviour is that t he more approximate t reatm ent of KLM 
and higher lines in lOelgoetz fc Pradhanl (|200ll ) appears to 
overestima te their importance, especia lly at higher temper- 
atures (see lOelgoetz fc Pr adhan 200l|, Fig. 3). This ovcrcs- 
timation leads to a cancelling effect, whereby the KLM and 
higher lines in the denominator of GD cancel out the effect 
of the KLL satellite lines in the numerator. Additionally, 
the present calculations allow the KLM and higher satellite 
lines to be included within the energy range they actually 
fall, which is in the redward section (i.e. the numerator of 
GD — with the x, y, and z lines) for some of the higher satel- 
lite lines. Thus, the satellite lines in these new calculations 
have an impact on the line ratio GD at temperatures well 
above the temperature of maximum abundance for the He- 
like ionisation stage. One practical consequence of this last 
statement is that essentially any spectral analysis of the He- 
like Ka lines requires the satellite lines to be treated in a de- 
tailed manner (unless the measured spectra are sufficiently 
well resolved so that the satellite lines can be readily dis- 
tinguished). Due to the level of detail and improved atomic 
data included in the present calculations, they are expected 
to be a significant improvement over previous work. 

While there are differences between the G ratios, as 
well as the GD ratios, predicted by each of the data sets, 
these differences are all less than 15 percent, which is 
within the typical 10-2 percent uncerta i nty reported for 
the R-matrix data (e .g. INahar et al"1l200ll ; |Pradhanlll985al ; 
IWhiteford et al.ll200lh . In order to understand these differ- 
ences, spectra were examined for a wide range of temper- 
atures. In general, spectra for all the elements and models 
considered were found to be in excellent agreement with 
each other, even when comparing results obtained from RM 
and DW data sets. The differences were all less than 12 
percent for strong lines, which include the w, x, y, and z 
lines, as well as most of the satellite lines. There were larger 
differences (up to ~50 percent) for some weak but barely 
visible satellite lines (like c), and even larger differences (up 
to ~150 percent) for some weaker satellite lines that do not 
contribute in any appreciable manner to the spectra. These 
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Figure 1. Top panel: GD and G ratios as a function of electron temperature for the Ca:DW and Ca:RM data sets. The temperature 
of maximum abundance (T m ~ 1.5 X 10 7 K) has been indicated with an arrow. Middle panel: Ratios of G and GD ratios computed 
from the various data sets. Bottom panel: Ratios of GD to G for each of the data sets. Insert on bottom panel: The same ratios of 
GD to G as the bottom panel, but plotted with a linear scale on the y axis to better resolve the behaviour. 



larger differences have very little impact on the spectra or 
the line ratios as the corresponding lines are quite weak. 

Two sample spectra for Fe, for which the disagreement 
in the ratios was among the largest, are presented in Figs. 2] 
and [5] As illustrated in the upper panel of Fig.|4j at an elec- 
tron temperature of 10 7 K, the overall agreement between 



the spectra computed with the various models is excellent. 
The data in the bottom panel of Fig. [4] indicate more pre- 
cisely where the largest discrepancies occur. One observes 
that the use of R-matrix data results in an increase of the 
z line and a decrease in the x line relative to the distorted- 
wave model. Additionally, the Fe:RM data set predicts a 
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Figure 2. Top panel: GD and G ratios as a function of electron temperature for the Fe:DW, Fe:RM. Fc:RM2, and Fe:RM+ data sets. 
The temperature of maximum abundance (T m ~ 4 X 10 7 K) has been indicated with an arrow. Middle panel: Ratios of G and GD 
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decrease in the y line, and an increase in the w line relative 
to the distorted-wave model; the Fe:RM2 data set predicts 
the same changes, but to a lesser extent. From this inspec- 
tion one can conclude that the agreement between the G 
ratios calculated from the Fe:DW and Fe:RM2 data sets is 



fortuitous because of a cancellation in the quantities that 
comprise the numerator and denominator of that ratio. On 
the other hand, the decrease in the x and y lines predicted 
by the Fe:RM versus the Fe:DW data set are larger than 
the corresponding increase in the z line. This overall reduc- 
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tion in the numerator of the G ratio, when coupled with 
the increase in the w line between the Fe:RM and Fe:DW 
data sets, results in the reduced G ratio calculated from the 
Fe:RM model at low temperatures. 

Fig. [5] which displays spectra at a much higher electron 
temperature of 4 x 10 8 K (which is approximately ten times 



higher than the temperature of maximum abundance for He- 
like Fe) again shows excellent agreement. An analysis of the 
bottom panel of Fig. [S] shows that both R-matrix data sets 
predict higher x and y lines, and a decreased w line, relative 
to the distorted-wave results. The net result of these differ- 
ences is the increased G and GD ratios displayed in Fig. [21 
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Separate calculations (not shown) indicate that the increase 
in the x line is due to slightly higher R-matrix recombina- 
tion rates rather than to sensitivity to the electron-impact 
excitation rates. This populating mechanism for the x line is 
consistent with the typical viewpoint in the literature (e.g. 
|Pradhanll985bT ). The y line, on the other hand, is sensitive to 
both electron-impact excitation and recombination rates at 
this high temperature; for this case the recombination rates 
are dominant in determining the population of the excited 
state, but the excitation rate is non-negligible as y is an in- 
tercombination line. While this temperat ure (T e = 4.0 x 10 8 
K) is above the peak of the DR hump (see lNahar et al ] |200ll . 
Fig. 5), it is still in a range where the resonances of the 
R-matrix cross section are important to the recombination 
rate. The high-temperature differences observed for the Ni 
G and GD ratios (Fig. [3}, for which only the recombina- 
tion rates were changed among the various models, have a 
similar explanation. Additionally, separate calculations (not 
shown) indicate that the differences in the w line are pri- 
marily due to differences in the electron-impact excitation 
data. The importance of excitation over recombination as a 
populating mechanism of the w line is expected since this 
transition is dipole allowed (e.g. |Pradhanlll985bl ). The net 
effect of these differences is the increase in the R-matrix G 
and GD ratios which is observed above the temperature of 
maximum abundance in Fig. [2] 

Despite the subtle differences in the spectra presented 
above, we emphasise that the discrepancies in the important 
lines are well within the uncertainties (20 percent) usually 
cited for R-matrix data. The disagreement in these spectra 
were among the largest seen in this study, which speaks to 
the excellent overall agreement between the RM and DW 
models. 

Lastly, it should be noted that the line positions for 
the KLM and higher satellite lines are a significant source 
of uncertainty in these calculations. While the accuracy of 
the line positions is estimated to be ~0.1 percent, a shift of 
that size could impact the spectra significantly by causing 
some of the strong KLM lines, which blend with the w line 
in this present work, to move sufficiently far such that they 
should be considered with the bulk of the KLL lines in the 
numerator of GD. This fact is underscored by the appear- 
ance of KLM and higher lines ~7 eV blueward of the w line 
in Fig. [4] when they should instead converge upon the w 
line. If some of these higher lying satellite lines do in fact 
blend with the x line, the impact would be a corresponding 
increase in the GD ratio. 



5 CONCLUSIONS 

New, more detailed calculations of the emission spectra of 
the He-like Ka complex of calcium, iron, and nickel have 
been carried out using atomic data from both distorted-wave 
and R-matrix calculations. Spectra from these calculations 
are in excellent agreement, and demonstrate that satellite 
lines are important to both the spectra and the GD ratio 
across a wide temperature range that includes temperatures 
significantly above the temperature of maximum abundance 
for the He-like ionisation stage. A major conclusion of this 
work is the need to include satellite lines in the diagnosis of 
He- like Ka spectra of iron peak elements in low density, col- 



lisional (coronal) plasmas, even at temperatures well above 
the temperature of maximum abundance. When the satel- 
lite lines are appropriately taken into account, the GD ratio 
remains an excellent potential temperature diagnostic. 

Another important application of the results presented 
herein, is in the well-known application of the G or GD ra- 
tio to ascertain the ionization state of a plasma. As shown 
in Figs. [1}J3] the GD ratio is far more sensitive to the ion- 
ization state at T < T m than the G ratio, by as much as 
a factor of 100. Therefore, it is imperative to calculate the 
GD values as precisely as possible at temperatures where 
the dielectronic satellite intensities are rapidly varying. Such 
conditions are known to occur in plas mas which are no t 
in coronal equilibrium, as discus sed by iPradhanl (|l985bh ; 
lOelgoetz fc Pradhanl (|200ll . l2004f ). Furthermore, it should 
be noted that, while this work does not consider the effect 
of satellite lines on the density sensitive diagnostic ratio R 
(R — (I(x) + I (y)) 1 1 (z)) , the effect of these lines is signifi- 
cant enough that they would need to be taken into account 
under conditions where R is used. This inclusion is war- 
ranted due to the manifestation of the satellites embedded 
within the Ka complex, and in many cases blended with the 
principal lines x, y and z. 

The excellent agreement between the spectra produced 
from R-matrix and distorted-wave data used in the models 
presented in this work bolsters confidence in both data sets. 
Any disagreement between the two sets of spectra would 
have indicated an error in the fundamental atomic data be- 
cause the IPIR approach has been shown to give good agree- 
ment with close-coupling approac hes when producing the 
fundamental rate coeff i cients (e.g. iBates fc Dalgarnol Il962l ; 
I Gabriel fc Pagetl 1 19721 ; ICowanl Il98d : iBadnell et al.l fl993h . 
The present work provides a more stringent test of this as- 
sumption by including those rate coefficients in a fully in- 
tegrated spectral calculation that takes into account several 
ion stages and includes the coupling between all of the im- 
portant atomic processes. 

The good agreement observed in this work reaffirms the 
fact that in highly charged systems, models based on data 
calculated from computationally less expensive distorted- 
wave methods can reproduce the results of models based 
on R-matrix data if the effect of resonances are taken into 
account as independent processes. This behaviour, however, 
is not expected to remain true for all conditions, especially 
when near neutral systems are prevalent. 

The results presented in this paper should be applica- 
ble to high-energy and high-resolution X-ray spectroscopy 
of laboratory and astrophysical plasmas. Astrophysical ob- 
servations of the Ka complex of high-Z ions, particularly the 
6.6-6.7 keV range of the Fe Ka, were expected to be made 
by the high-resolution X-ray satellite Suzaku, but could not 
be performed due to instrument failure. It is, however, ex- 
pected that these calculated results would be valuable in 
future X-rays missions such as the recently planned joint 
ESA-NASA International X-ray Observatory. 
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Figure 5. Top panel: Calculated Fe spectra at T e = 4.0 X 10 8 K, N e = 10 10 cm -3 for each of the Fe data sets. Bottom panel: 
Functions of the spectra presented for comparison. Again there is excellent agreement between all spectra. The largest uncertainties are 
in very weak satellite lines, as well as the w and x lines. 
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